The prediction method of similar cycles 
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ABSTRACT 

The concept of degree of similarity (ry) is proposed to quantitatively describe the similarity 
of a parameter (e.g. the maximum amplitude i?max) of a solar cycle relative to a referenced one, 
and the prediction method of similar cycles is further developed. For two parameters, the solar 
minimum (i?min) and rising rate (/?a), which can be directly measured a few months after the 
minimum, a synthesis degree of similarity {rjs) is defined as the weighted- average of the rj values 
around i?min and /3a with the weights given by the coefficients of determination of i?max with 
i?min and /3a, respectively. The monthly values of the whole referenced cycle can be predicted 
by averaging the corresponding values in the most similar cycles with the weights given by the 
77s values. Cycles 14 and 10 are found to be the two most similar cycles of Cycle 24. As an 
application. Cycle 24 is predicted to peak around January 2013 ± 8 (months) with a size of 
about i?max = 84 ± 17 and to end around September 2019. 



Subject headings: Sun: activity — Sun: general — sunspots 



1. Introduction 

Solar activity prediction is important for both 
space weather and solar physics. Since solar activ- 
ity is the driver of various phenomena in the near- 
Earth environment, knowing the future level of so- 
lar activity in advance can reduce some of the haz- 
ards for high-tech equipment on which our mod- 
ern society dependents. Successful predictions 
could p rovide some constraint s on sol ar dynamo 
models (ICameron fc Schusslei|[20 08: P esnellll2008t 
Wang et al.ll2009at ICuo et al.ll2010l: lDull2011al) . 

Various techniques have been used in the past 
to predict solar activity, of which some were purely 
statistical and others were related to physics 
( Kand 120071: IZhan-leDu et all 120081: |Pesnellll2008t 



netic field that gener ates the tor oidal field for the 
next cycle ((Schatten et al.lll978l ). Solar dynamo 
models have recently been attempted, but they 
are as yet im mature to predict the associated 



solar activity (iDikpati et al, 
20061: lPesnellll2008l) . 



2006t iTobias et al.1 



A prominent feature in solar activity is the 
so-called Waldme ier effect, where stronger cycles 
tend to rise faster ( Waldmeieill 1 9391 : Hathawav et al 



Messerotti et al.l 120091) . Geomagnetic precursor 
methods have attracted more attention due to 



successes in Cyc l es 20-22 (IOhi|ll966l: iKand 12007 
Du et all l2009al : IDu fc Wane! boilh . They are 
based on a solar dynamo concept that the geo- 
magnetic activity during the declining phase of 
the preceding cycle or at the minimum provides 
an approximate measure of the poloidal solar mag- 



2002|:IDu et al.ll2009bl ). This effect implies that the 
magnetic energy in a solar cycle has a tendency of 
stability and that str onger cycl es need less time to 
release their energies (lDull2006l) . Correspondingly, 
the maximum amplitude (i?max) of a solar cycle 
is well correlated with the growth rate of activ ity 
in the early phase ( Cameron fc Schiissler 20081 ). 

As a new solar cycle is ongoing, its shape can 
be well described by simple functions contain- 
ing a few parameters jElling fc Schwen tck 199^ 
Hathawav et~al]|l994 lLilll999l: IVolobuevi,2009i: .Dul 



2011br ) . Before the arrival of the peak of a solar cy- 
cle, only two parameters are known: the preceding 
minimum (i?min) and the rising rate (/3a), defined 



1 



as the ratio of the increment of the activity (i?z) 
and the elapsed time duration. The two parame- 
ters can be taken as indicators for the subsequent 
amphtude (i?max)- 

Solar cycles that have approximatively the 
same i?max tend to have similar shapes, and 
these cycles are therefore cal l ed 's i milar cycles' 
dWaldme ier' 'l936|- iGleissberd Il971t IWang et al 



20oi) . iGlcissbcr j ([1971^ used the concept that 



similar cycles tend to ha ve similar cycle l engths 
and similar decline times ( Gleissberg 19731) to es- 
timate the epochs of the start and end of Cycle 21 
by averaging the cycle length s and decline t i mes o f 
similar cycles, respectively. IWang fc HaJ (|l997h 
developed the co ncept of similar cyc les and used it 
to predict i?rna.x (IWang et al.ll2009bl) and monthl y 
values of |wang et al.ll2002t iMiao et al.ll2008[ ). 
The predicted value of Rz for month i is taken as 
the average of the corresponding values of similar 
cycles for the same month from the start of the 
cycles, and the standard deviation of these val- 
ues is taken as the prediction error for the same 
month. 

The data and parameters used in this study are 
shown in Section [21 Some of the correlations be- 
tween these parameters are analyzed in Section 
[3l and in Section |4l we employ the two parame- 
ters of i?inin and /3a to further develop the concept 
of a similar cycle and its application in predic- 
tion. In Section 14.11 a quantity, the "degree of 
similarity is proposed to describe the sim- 

ilarity of a parameter in past cycles relative to 
a referenced one (the predictor for a cycle to be 
predicted). For two predictors (i?min, Pa), in Sec- 
tion [421 a synthesis degree of similarity (77s) is de- 
fined as the weighted-average of the corresponding 
77 values (77R,, rjp) with the weights given by the 
coefficients of determination of -Rmax with i?min 
and Pa (j'r, rj), respectively. Then, the i?niax 
value for the current cycle (24) can be predicted 
as the weighted-average of the -Rmax values of the 
five most similar cycles with the weights given by 
rjs ■ The same technique is used in Section 14.31 for 
each month from the start of the similar cycles 
to obtain the monthly values (shape) of Cycle 24. 
The predictive power of the similar-cycle predic- 
tion method is tested on different months from the 
start of Cycle 24 in Section 14.41 The results are 
briefly discussed and summarized in Section [5l 



2. Data and cycle parameters 

The data used in this study are the smoothed 
monthly mean Ziirich relative sunspot numberf0 
(i?z) of the more reliable data since Cycle 8 (up 
to November 2010). Some parameters of the solar 
cycle are listed in Table [1] where i?max and i?rnin 
are the maximum and minimum amplitudes of a 
solar cycle, respectively. Here /3a is the rising rate 
at Am months entering into the cycle. 



/3a = Ai?z/Am, 



(1) 



where AR^ — Rz{Am) — i?min is the increment of 
Rz from i?inin in the time interval of Am (months). 
At the present time. Am = 24 (months) for the 
data available in the current cycle (24). S is 
the skewness of the cycle, a measure of symme- 
try about its maximum. 



(/V- 1)^3 



S = 



(2) 



where y = i?z is the data series for a solar cy- 
cle, y is the mean, N the number of data points, 
and a the standard deviation. Positive (negative) 
skewness indicates that the distribution is skewed 
to the right (left), with a longer tail to the right 
(left) of the maximum. 
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Fig. 1. — (a) /3a as a function of Am for Cycles 8- 
24. The dashed lines indicate /3a for Cycles 14, 10, 
17, 16 and 12. The thick solid line indicates the 
averages for the data available, (b) Correlation 
coefficient (r) between i?max and /3a. 

In Table [1] x and a are the average and stan- 
dard deviation of the parameters for Cycles n = 8- 
23, respectively, r is the correlation coefficient of 
the parameter with time (temporal variation), and 
CL is the confidence level. 

It is shown in Tabled that i?min7 /3a and i?max 
all show an increasing trend with time (r > 0) at 



^ http: / / www.ngdc.noaa.gov/stp / spaceweathcr.html 
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Table 1: Parameters since cycle n = 8. 
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7.3 


2.85 


3.31 
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0.28 
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10.6 


1.10 


1.25 


131.9 


0.51 


10 


3.2 


1.24 


1.37 


98.0 


0.14 


11 


5.2 


2.48 


2.63 


140.3 


0.52 


12 


2.2 


1.65 


1.73 


74.6 


0.09 


13 


5.0 


2.32 


2.38 


87.9 


0.35 


14 


2.7 


1.28 


1.37 


64.2 


-0.09 


15 


1.5 


1.94 


1.95 


105.4 


0.33 


16 


5.6 


1.29 


1.73 


78.1 


0.02 


17 


3.5 


1.46 


1.79 


119.2 


0.12 


18 


7.7 


2.11 


2.47 


151.8 


0.16 


19 


3.4 


4.07 


4.80 


201.3 


0.30 


20 


9.6 


1.94 


2.42 


110.6 


0.03 


21 


12.2 


2.32 
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164.5 


0.13 


22 


12.3 


3.88 


4.54 


158.5 


0.16 


23 
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2.14 
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0.30 
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1.7 


0.75 


1.03 
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Standard deviation. 

Correlation coefficient of temporal trend. 
The confidence level (%). 
At the current state. 



the confidence level (CL) around 80%, and that 
S used to be positive {S = 0.21), meaning that 
the decline time tends to be longer than the rise 
time for most solar cycles. However, this asym- 
metry seems to decrease as can be seen from the 
decreasing temporal trend (r = —0.32) in S. 

The /?a values for Cycles 8-24 are shown in 
Fig. [IJa) . It is seen in this figure that /3a varies ap- 
proximately linearly with Am at the early phase of 
the cycle (about Am < 26) and has decreases since 
then. The correlation coefhcient (r) between i?max 
and /3a varies with an increasing trend (Fig. [T]d) 
and r > 0.75 when Am > 21. First, we take 
Am = 21 (months) as an example to describe the 
method of similar cycles. Then, this method is 
applied to Am = 18, 19, • • • , 24 (months) for pre- 
dicting Cycle 24. 



3. Correlation betvifeen some parameters 

Of the parameters describing the solar cycle, 
only Rmin and /3a are known before the timing 
of i?max- Therefore, both i?min and /3a can be 
taken as indicators for the subsequent i?max- Fig- 
ure[2][a) shows the scatter plot of -Rmax against 
-Rmin with the least-square-fit linear regression 
equation given by 



Rn 



91.3 + 4.94i?min,o- 33.1, 



(3) 



where a is the standard deviation. It is well known 
that i?max is weakly correlated with the preceding 
^min (''R = 0.47), at the 92.9% confidence level, so 
that a small i?min tends to b e fo llowe d by a weak 
cycle (|Hathawav et al.ll2002l: 111 120091) . However, 
since this correlation is weak, one can hardly ob- 
tain an accurate pr ediction of i?max by i?min alone 
ton fc Wandl2010[) . 




Fig. 2. — Scatter plots of (a) i?,nax against i?min, 
(b) -Rmax against ^a, (c) S against i?max, (d) 5' 
against /3a. 

Figure[2Ib) shows the scatter plot of -Rmax 
against /3a. The best linear relationship between 
them is 

-Rmax = 54.4 + 31.98/3a, a = 24.8. (4) 

One can see that -Rmax is well correlated with 
the rising rate /3a {rp = 0.75) at the 99.9% con- 
fidence level. Therefore, /3a is a good predictor 
for the subsequent -Rmax- Substituting the value 
of /3a (0.75) for Am = 21 into this equation, 
the peak size of Cycle 24 can be estimated to be 
i?,„^^(24) = 78.3 ± 24.8 (triangle). 
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The scatter plots of S against Rmux and S 
against /3a are shown in Figs, ^c) and (d), re- 
spectively. The correlation coefficients involved 
are listed in Tabled 

Table 2: Cross-correlation Coefficients (r). 
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Rmax 
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0.47 


92.9% 
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0.75 
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S 


-^max 


0.41 


88.5% 


S 


/3a 


0.24 


62.0% 


s 


^min 


0.10 





4. Similar cycles 

It should be pointed out in Fig. [5Jc) that S 
is positively correlated with i?max (»" = 0.41), so 
that similar cycles (with approximatively the same 
J?max) tend to have similar shapes, which is a key 
point in the concept of similar cycles although r 
is not high. Because if S were uncorrelated with 
^^max, the concept of similar cycles would not have 
been used any longer. However, as i?max is un- 
known in advance, the shape of an upcoming cycle 
cannot be described directly from the measure- 
ment of i?max, which is usually estimated by an 
appropriation technique. 

Because i?max has a much higher correlation 
coefficient with /3a (r = 0.75) than with i?min 
(r = 0.47) and the correlation coefficient of S with 
/?a {r — 0.24) is slightly higher than that of S with 
Rmin (r = 0.10), /3a is & much better predictor for 
.Rmax than i?niin- Howcvcr, the two predictors of 
/3a and i?inin are expected to provide useful infor- 
mation (size and asymmetry) on the subsequent 
cycle, because both of them can be directly mea- 
sured a few months after the minimum. 

4.1. Degree of similarity 

For a parameter x — {xi, i — 1, 2, • • • , m}, the 
probability density of the normalized deviation. 



Aa 



(5) 



approximately satisfies a normal distribution 
(Fig.ISi), 

p(E,) - ^e-^'/2, (6) 



/27r 



where x and ctx are the mean and standard devia- 
tion, respectively (Table [Ij . The probability that 
Axi falls in the range [Azr, Axn] is 

AP(n) ^ e-*'/"di, (7) 



27r 

where referenced value (a predictor for Cy- 

cle n-c — 24) and Xn the value in the past cycle 
(n). 
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Fig. 3. — (a) Gaussian distribution, (b) i?min 
(soHd) and t^r (dashed), (c) /3a (solid) and -qp 
(dashed). The numbers indicate the similar cycles 
and their orders in brackets. 

A small value of |AP(n)| indicates that Xn{= 
X + Axn) is close to Xr{= x + Ax,-), in which case 
Cycle n is called a similar cycle of around x^- 
The smaller the |AP(n)| is, the more similar the 
two parameters are. Because < |AP(n)| < 1 for 
either Xr > Xn or x^ < Xn, we define 1 — |AP(n)| 
as a measure to describe the "degree of similarity" 
around x,-. 



ri{n) 



1 



1 



/2n 



(8) 



If ri{n) = 1, the two parameters are identical, 
a^n = Xi- (100% similarity); if rj^n) = 0, there is 
no similarity between them. In application, some 
(five in this study) of the largest values of 77(71) 
are selected and their cycles are called the similar 
cycles of around x,. 

4.2. Similar cycles used in predicting Pmax 

Figure|31^b) shows the values of Pmin (solid) for 
Cycles n = 8-23, the referenced one (Pmin(24) = 
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1.7, triangle) for Cycle n,- — 24, and the val- 
ues of r] (?7R, dashed) around i?min(24). The 
five biggest values of t^r (asterisks) are in Cycles 
TT-R — 15,12,14,10 and 19 (in that order), which 
are called the similar cycles of Cycle 24 around 

Figure 3(c) shows the values of /3a (sohd), the 
referenced one (/3a(24) — 0.75, triangle), and the 
values of r] {r^jj, dashed) around /3a (24). The 
five biggest values of 77^ (asterisks) are in Cycles 
np = 9,10,14,16 and 17 (in that order), which 
are called the similar cycles of Cycle 24 around 
/3a. These similar cycles (n^) are partly different 
from the above ones around i?niin («r). Since the 
correlation coefficient of i?,nax with /3a {rp = 0.75) 
is much higher than that (rR — 0.47) of i?,nax with 
-Rmini f^p should bc morc reliable than tt-r for de- 
scribing the information of Cycle n,.. 

From 77R and 77^ , a synthesis degree of similarity 
around both i?niin and /3a can be defined as 



77s (n) 



(9) 



'/3 



here the weights are taken as the coefficients of de- 
termination, — and — for i?niin and 
/3a, respectively. The reason for such a selection is 
that about (r|) of the variations in i?max can 
be explained by the correlation of i?max with i^min 

(/3a). 
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Fig. 4. — (a) i?max (sohd) and 77s (dashed), 
(b) The values of i?z for similar cycles Us = 
14, 10, 17, 16 and 12. (c) The predicted monthly 
values of Rz for Cycle 24 with error bars. 

FigureSfa) shows the values of i?niax (solid) and 



77s (dashed). The five biggest values of 77s (aster- 
isks) are in Cycles 7is = 14, 10, 17, 16 and 12 (in 
that order), which are called the similar cycles of 
Cycle 24 around both i?i„in and /3a. From these 
values, the peak size of Cycle 24 can be predicted 
as the weighted-average of those of the similar cy- 
cles. 



-^max(24) 

W{i) 



ELl^max(n.(i))W(») 



(10) 



where the weights W{i) are taken as the (synthe- 
sis) degrees of similarity 77s(7is(i)). The more simi- 
lar a cycle is to 7ir, the more its weight it should be 
included. The standard deviation is correspond- 
ingly defined as 



.(24) 



'Etl[^max(ns(«)) - i?„,ax(24)]2W^(l) 



(11) 

According to the above equations, the peak size of 
Cycle 24 is predicted to be i?max(24) = 86.8 ± 19.6 
(triangle). 

4.3. Similar cycles used in predicting 
monthly values 

Now, we use the above technique to predict the 
monthly values (shape) of Cycle 24. Figure|Ub) 
shows the monthly values of for the similar cy- 
cles {us) from the starting points of the cycles. 
The i?z value of the t*'' month in Cycle 24 is pre- 
dicted as the weighted-average of the correspond- 
ing Rz values for the same month in the similar 
cycles. 



Rz{t) = 



E.Li^z(t,r7s(7))m*) 



Its standard deviation is 



(12) 



O-z(t) 



'j:limt,ns{t))-Rz{tWW{z) 



(13) 

The results since November 2008 are shown 
in Fig. m^c). From these results, we can obtain 
the maximum value (83.0), the standard deviation 
(16.7), the rise time (Ta — 58 months), and the cy- 
cle length (130 months). Cycle 24 is therefore pre- 
dicted to peak around September 2013 with a size 
of about i?,nax — 83.0 ± 16.7, and to end around 
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September 2019. This size is near that (86.8) in 
section 221 and that (78.3) in Section |3l It should 
be pointed out in Fig.SJc) that the shape near the 
peak is rather flat. Besides the highest peak (83.0) 
in September 2013, there are two weak peaks pre- 
ceding it: the first is around June 2012 (81.1) and 
the second is around January 2013 (80.1). If this 
result is true, it implies that Cycle 24 will have 
multiple peaks with the last peak being higher. 
On average. Cycle 24 may probably peak around 
January 2013 ±8 (months). 



4.4. Prediction results at 
months 



Am 



18-24 



Since /3a is a temporal variable of Am, the result 
derived from /3a may depend on Am, similar to 
the prediction by a simp le function to describe the 
shape of the solar cycle (jHathawav et al ] |1994I: Id^ 
2011br ). In this section, we apply the above tech- 
nique to the current state (Am — 24), as shown in 
Fig-El The similar cycles at Am = 24 (rig =14, 10, 
12, 17 and 15) are slightly different from those at 
Am = 21 (ns =14, 10, 17, 16 and 12) in Fig.lH In 
Fig.jSJa), the peak size of Cycle 24 is predicted to 
be i?,„ax(24) = 91.6 ± 20.1 (triangle) based on the 
-Rmax values of the similar cycles, slightly higher 
than that (i?max(24) = 86.8 ± 19.6) in Fig.HJa). 
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Fig. 5. — Similar to Fig. [H but using the /3a val- 
ues at Am = 24 (month), (a) i^max (solid) and 
77s (dashed), (b) The values of i?z for similar cy- 
cles 71s = 14, 10, 12, 17 and 15. (c) The predicted 
monthly values of i?z for Cycle 24 with error bars. 

In Fig. EKc), the highest peak is i?max = 86.5 ± 
21.4 in January 2013 (from the predicted monthly 



Rz in Cycle 24), which is near the second peak 
in Fig. Sfc). In addition, there are two shoulders: 
one in July 2012 (83.8) and another in Septem- 
ber 2013 (84.5). The former is near the first peak 
(June 2012) and the latter is just the third peak 
(September 2013) in Fig. (He). 

Using the same technique, we tested the predic- 
tive power of this method at Am = 18-24 month. 
The results are listed in Table [31 where the third 
column indicates the similar cycles for a given Am 
(first column); the fourth- sixth columns corre- 
spond to the first -third peaks, respectively; and 
the last column is the ending time of Cycle 24. The 
last row shows the relevant averages. It is seen in 
this table that there are not significant differences 
in the results as the cycle progresses although the 
similar cycles may have a small difference (us) for 
different Am. The three peaks for different Am 
are near each other either in size or in date. If 
the maximum average size and the middle date 
is taken as those for the peak of the cycle. Cycle 
24 is predicted to peak around January 2013 ± 8 
(month) with a size of about i?max = 84 ± 17. 

5. Discussions and Conclusions 

A concept, the degree of similarity (rj), is pro- 
posed to quantitatively describe the similarity of a 
parameter of a solar cycle relative to a referenced 
one, and the prediction method of similar cycles 
is further developed. The degrees of similarity are 
used as the weights in the weighted-average of the 
values of a parameter in the similar cycles so as 
to obtain a predicted one in the referenced cy- 
cle, where we have considered the fact that more 
weights should be paid to the cycles that are more 
similar to the referenced one. 

In this study, we used two predictors, the pre- 
ceding minimum (Rmin) and rising rate (/3a), to 
define a synthesis degree of similarity (rjs) by aver- 
aging the corresponding rj values with the weights 
given by the coefficients of determination of i?max 
with i?„iin and /3a, respectively. From this method. 
Cycle 24 is predicted to peak around January 2013 
± 8 (months) with a size of about i?inax = 84± 17 
and to end around September 2019. This result 
is s lightly lower than th at (100.2 ± 7.5) predicted 
by IWang et"aL ( 2009b[ ) based o n their similar - 
cycle method. It is near 80 ± 30 (ISchatte3 l2005l). 
~ 80 (|Kitiashvih fc Kosovich"evll2008l ). and 85 
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Table 3: Prediction Results for Am = 18, 19, • • • , 24 (months). 



Am 


/3a 




-Rmax(lst peak) 


i?niax(2nd peak) 


-Rmax(3rd peak) 


ending 


18 


0.77 


10,17,14,16,15 


85.9±20.5(6/2012) 


89.4±19. 2(12/2012) 


86.4±15.8(8/2013) 


2/2019 


19 


0.77 


10,17,14,16,12 


81. 8±22. 7(6/2012) 


80.7±20.6(1/2013) 


83.5±16.8(9/2013) 


9/2019 


20 


0.75 


10,14,17,16,12 


81.4±22. 6(6/2012) 


80.4±20.6(1/2013) 


83.2±16.7(9/2013) 


9/2019 


21 


0.75 


14,10,17,16,12 


81. 1±22. 6(6/2012) 


80.1±20.5(1/2013) 


83.0±16.7(9/2013) 


9/2019 


22 


0.81 


14,10,17,12,16 


81. 0±22. 5(6/2012) 


80.0±20.5(1/2013) 


82.9±16.7(9/2013) 


9/2019 


23 


0.93 


14,10,12,17,16 


80. 9±22. 5(6/2012) 


79.9±20.5(1/2013) 


82.8±16.7(9/2013) 


9/2019 


24 


1.03 


14,10,12,17,15 


83. 6±22. 6(7/2012) 


86.5±21.4(1/2013) 


84.5±16.5(9/2013) 


8/2019 


X 






82. 2±22. 3(6/2012) 


82.4±20.5(1/2013) 


83.8±16.6(9/2013) 


9/2019 



(jjiang et al.l 120071 ) based on polar field or solar 
dynamo models, but much lower than ~ 167 



(jPikpati et al.l 120061 ) based on a modified flux- 



transport dynamo model. 

It should be pointed out from Section 
that (i) considering a single parameter i?min, the 
most similar cycle (to Cycle 24) is n = 15 since 
i?min(15) = 1.5 is close to i?min(24) = 1.7; (ii) 
considering another single parameter /3a, the most 
similar cycle is n = 9 as /?a(9) = 1.10 is close to 
/3a(24) — 0.75; and (iii) considering both parame- 
ters i?inin and /3a, thc most similar cycle is n = 14. 
In this study, the similar cycles are selected from 
two parameters (i?min, /3a) to avoid an acciden- 
tal error caused by a single parameter. In terms 
of the synthesis degree of similarity (77s), Cycle 
15 is not a similar cycle since /3a (15) = 1.94 is 
not so near to /3a(24), although i?min(15) is close 
to i?min(24), and Cycle 9 is not a similar cycle 
as i?min(9) = 10.6 is far from i?inin(24) although 
/3a(9) is close to /3a(24). Therefore, the most simi- 
lar cycle (n = 14) is that with maximum rjs which 
is the synthesis effect of both i?i„m(14) = 2.7 and 
/3a^(14) = 1.28, although Cycle 14 is not the most 
similar cycle with either Rmin or /3a alone. Be- 
cause /3a varies with the progression (Am) of the 
cycle, the results from both i?,nin and /?a may have 
a slight differences (Table [3]). 

Near the time of a solar minimum, geomagnetic 
activity is a much better indicator for the ensu- 
ing maximum amplitude (-Rmax) of the solar cy- 
cle ( Ohll 1966 ) than the solar minimum (i?rnin)- 
This study shows that the rising rate (/3a) at the 
early phase of a solar cycle is also a good indica- 
tor for the subsequent i?max- This parameter has 
the advantage that it only needs the series it- 
self. It reflects the initial physical process of solar 



magnetic activities. 

Similar cycles are usually defined as those 
whose values of i?max (or other parameters) satis_ 
a given condition dGleissberg 1971; ..Wang fc .Ha 
19971 : IWang et al.ll2002l i2009bt iMiao et al.ll2008l ). 

)l < A, (14) 

where i?max("-r) is the referenced value and A is 
a given limit {e.g., 10 or 20). By averaging the 
values of a parameter (rise time, decline time or 
monthly i?z and so on) in these cycles, the cor- 
responding one in the predicted cycle could be 
obtained. This technique has considered the cy- 
cles that have similar i?niax to i?max("-i)- However, 
how close these cycles are to the referenced one has 
not been considered, as a simple averaging method 
{W{i) = 1) was used in this technique. Besides, 
^max(f^r) IS usually unkuowu, so it also needs to 
be predicted. The error in i?max("T) may probably 
propagate into the other parameters that are de- 
rived from it. In contrast, this study used the two 
directly measured parameters (i?min, /3a) to derive 
all the information of a cycle to be predicted (24). 

The central idea of a similar cycle is that solar 
cycles with approximatively the same i?max tend 
to have similar shapes, which is empirical and has 
not been studied physically as far as we know. It is 
shown in Fig.[2Jc) that S is weakly correlated with 
^max {t — 0.41), which is a key point in the con- 
cept of similar cycles. Under this relationship, we 
can use the concept that solar cycles with approx- 
imate i?max might have similar rise times, cycle 
lengths and cycle shapes, so we can estimate it in 
the referenced cycle by averaging the correspond- 
ing values in similar cycles. However, since i?max 
is unknown in advance, the shape of an upcom- 
ing cycle cannot be estimated directly from i?max- 
In this study, we used both the rising rate (/3a) 
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and solar minimum (i?rnin) to select the similar 
cycles, which is based on the correlations of .Rmax 
with both /3a {rp ~ 0.75) and i?min {tr = 0.47). 
If one uses other parameters (e.g., the preceding 
decline time) to find some similar cycles, the cor- 
relations between iimax and these parameters are 
also needed (even if they are not strong). 

One advantage of the similar-cycle method is 
that it does not involve the details of a phys- 
ical process. The actual physical process may 
be rather complicated; its dynamical mechanism 
is not very clear at present and may not be de- 
scribed by a simple linear or non-linear relation- 
ship. Whatever the process is, a similar process 
may likely occur if the levels of activity are sim- 
ilar, which is what the similar-cycle method can 
and wants to do. This is similar to treating the 
complex process as a piecewise function. 

The main points in this study can be summa- 
rized as follows. 

1. A concept, the degree of similarity (77), is 
proposed to quantitatively describe the sim- 
ilarity about a parameter for a solar cycle 
relative to a referenced one. 

2. For two parameters, the preceding mini- 
mum (i?min) and rising rate (/3a), a synthe- 
sis degree of similarity {%) is defined as the 
weighted-average of the corresponding ones 
with the weights given by the coefficients of 
determination of i?max with i?niin and /3a, re- 
spectively. 

3. The prediction method of similar cycles is 

further developed with the weights given by 
the (synthesis) degrees of similarity. 

4. As an application of this method. Cycle 24 is 
predicted to peak around January 2013 ± 8 

(months) with a size of about i?max = 84±17 
and to end around September 2019. 
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